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Aspects of the dynamical glass transition are considered 
within a mean field spin glass model. At the dynamical tran- 
sition the the system condenses in a state of lower entropy. 
The difference, the information entropy or complexity, is cal- 
culated by analysis of the metastable (TAP) states, ft in- 
creases for lower temperatures, showing that more and more 
free energy barriers cannot be overcome on accessible times 
scales. Near T c the total glassy free energy lies below the para- 
magnetic one, which makes the transition unavoidable from a 
thermodynamic point of view. The multitude of glassy states 
implies an extensive difference between the average specific 
heat and the derivative of the average internal energy. 
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The structural glass transition is said to occur at the 
temperature T g where the viscosity equals 10 Poise. 
The precise nature of this transition is not fully under- 
stood. Experimentally one often determines the entropy 
from the specific heat data by integrating C/T from a 
reference temperature in the liquid phase down to T. At 
zero temperature a residual entropy is found, in contra- 
diction with the third law of thermodynamics. 

Alternatively, a glass can be seen as a disordered 
solid. In this description the glass transition is a tran- 
sition to a state with much smaller entropy than the 
liquid. As the free energy then becomes larger, it is 
not immediately clear from thermodynamic considera- 
tions why the system gets captured in such a state. 
The answer is that the condensed system then has lost 
the entropy of selecting one out of the many equivalent 
states; this part of the entropy is called the informa- 
tion entropy or complexity X. m @] Its origin can be un- 
derstood as follows. When the relevant state a has a 
large degeneracy J\f a = exp(2" a ), the partition sum yields 
Z = Efc exp(-/3F b ) « M a exp(-/?F a ), so F = F a -TJ a is 
the full free energy of the system. The entropy loss arises 
when the system chooses the state to condense into, since 
from then on only that single state is observed. H As the 
total entropy S — S a + T a is continuous, so is the total 
free energy. If the weights p a of the states a are fixed at 
the transition, jj| one can assume that the free energy 
difference between the condensed phase and the liquid is 
positive and grows quadratically below T c . This explains 
the well known discontinuity in quantities such as the 
specific heat. However, it is a bit unsatisfactory that this 
higher free energy branch describes the physical state. 

Recent work on glasses mainly has been performed 



from a dynamic approach. Mode coupling equations 
starting from the liquid phase have been studied for some 
while. Q At a critical temperature T c > T g a dynamic 
phase transition has been reported. The presence of this 
transition has been questioned, however. j5j 

We wish to analyze these and related questions within 
a relatively well understood spin glass model, the mean 
field p-spin interaction spin glass, defined by the Hamilto- 
nian (|l|) . Kirkpatrick and Thirumalai |(| pointed out that 
for the case p = 3 there is a close analogy with models for 
the structural glass transition. Moreover, its properties 
are quite insensitive of the value of p as long as p > 2. 
The spherical limit of this model has received some atten- 
tion recently. Its static version was solved by Crisanti and 
Sommers. They find a static first order transition to 
a state with one step replica symmetry breaking (1RSB) 
at a temperature T g . The dynamics of this system was 
studied by Crisanti, Horner and Sommers (CHS) || and 
Cugliandolo and Kurchan ||. Both groups find a dy- 
namical transition at a temperature T c > T gi which can 
be interpreted on a quasi-static level as a 1RSB transi- 
tion. For T < T c one of the fluctuation modes is massless 
("marginal"), not unexpected for a glassy state. At T~ 
there is a drop in the specific heat. Integrating C/T to 
define the "dynamical" entropy yields a free energy that 
exceeds the paramagnetic one quadratically. Ref. pQ ] 
discusses the interpretation of metastable states ("TAP- 
states" ) in this system. The statistics of those states was 
considered by Crisanti and Sommers (CS). [jllj They re- 
produced the dynamical free energy obtained in || . As a 
side remark, let us draw attention to the model with both 
random pair couplings and random quartet couplings. It 
is an exactly solvable spin glass model with infinite order 
replica symmetry breaking, p2| which remains exactly 
solvable after quantization of the spherical spins, pjfl 

The present author recently proposed a replica calcula- 
tion as a short cut for the long time result of the dynami- 
cal analysis, [jlij He took the principle point of view that 
a dynamical transition will automatically get trapped in 
a state with diverging time scale, if present. In a replica 
analysis this led him to replace Parisi's stationarity cri- 
terion for the location of the breakpoint of the 1RSB so- 
lution by a marginality criterion. In this way the correct 
dynamical transition temperature and shape of the order 
parameter function were reproduced from replica's for a 
number of model systems. However, as the free energy is 
minimized in this procedure, it was found to lie near T c 
below the paramagnetic value and to have a larger slope. 



1 



Though this is exactly what one expects at a first order 
phase transition, such a replica result cannot be trusted 
without independent verification and interpretation. 
For a system with N spins we consider the Hamiltonian 
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with independent Gaussian random couplings, that have 
average zero and variance J 2 p\/2N P ~ 1 . The spins are 
subject to the spherical condition ^ - Sf = N. At zero 
field the replica calculation gives the free energy 



F 

N 



(3J 2 
4 



PJ 2 
4 



£q P 



(2) 



T 

2x, 



log(l - £q) 



2x 



log(l - q) 



where £ = 1 — x. The first term describes the paramag- 
netic free energy. Here and in the sequel, we omit the 
T = oo entropy. It is a constant, only fixed after quan- 
tizing the spherical model, (T^| that plays no role in the 
present discussion. For the marginal solution q is fixed by 
equating the lowest fluctuation eigenvalue to zero, which 
gives p(p - l)f3 2 J 2 q p ~ 2 (l - q) 2 /2 = 1. The condition 
dF/dq — then yields x = x{q) = (p — 2)(1 — q)/q. 
This dynamical transition sets in at temperature T c — 
J{p(j>— 2) p ~ 2 /2(p— ljP -1 } 1 / 2 where x comes below unity. 

The metastable states of this system, labeled by a = 
1, 2, • • • ,VV, have local magnetizations m° = (Si) a . At 
H = their free energy F a = FxApinT-i) is a minimum 
of the "TAP" free energy [H 
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where now q = (1/-/V) ^\ mf is the self-overlap. Palmer 
Q calls the metastable states "components". Accord- 
ing to the Gibbs weight their probability of occurrence is 
Pa = exp(-pF a (T))/Z, with Z = J2a e M-PF a ). Given 
the p^s one can calculate the "component averages" such 
as F = J2 a Pa F a and S = J2 a Pa^a- For a given sample 
they are the objects one finds when repeating the mea- 
surement very often, and averaging over the outcomes. 

The nice thing of the present model is that many 
questions can be answered directly. After setting 
dFxAp/drrii = 0, we can use this equation to express 
F a in terms of q a alone. This gives F a = Nf(q a ) where 



/(?) 



4 
Tq 



l + (p-l) q P-(p-2)qP- 1 } 



T 

- log(l - q) 



(4) 



Since F a only depends on the selfoverlap q a , it is self- 
averaging. In the paramagnet one has m., = q — 0, so 



both eqs. (||) and (||) reproduce the replica free energy 
F/N = —(3 J 2 /4. From the replica analysis we know 
that at T~ the value of q is q c — (p — 2)/(p— 1). The 
component free energy F a — Nf(q c ) exceeds the free 
energy —f]J 2 N/A of the paramagnet. As expected from 
experimental knowledge on glasses, the internal energy 
is continuous. The difference is solely due to the lower 
entropy, S a = —N(3 2 J 2 /4 — I c , where 



J, =N[ ilog(p-l) + --l 
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is the value of complexity at the transition point. 

This discussion shows explicitly that the entropy of the 
component the system condenses in, as well as its average 
S, is much smaller that the entropy of the paramagnet. 
In the quantized system S will vanish at T = 0. fl2]| 

A simple calculation shows that the dynamical free en- 
ergy of CHS and CS, and the marginal replica free energy 
obtained from eq. (||) JlJ] have the following connection 
with the component average free energy F — Nf(q): 
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Since x = 1 at the transition, both are continuous there. 

Due to the difference between (Q) and (g) we have de- 
cided to redo the analysis of the TAP equations. Hereto 
we consider the generalized partition sum 
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For ii = we thus calculate the total number M of 
TAP-states, while for u = 1 we consider their par- 
tition sum. The sum over the TAP states can be 
calculated using standard approaches. (ll| One goes 
to integrals over all Wi-values of exp(— uf3FTAp{mi)) 
Yl^S (d/3F T Ap/dmi)det(d 2 PFTAp/dmidrrij) @ and in- 
troduces conjugate fields for writing the delta functions 
and the determinant as exponential integrals. After av- 
eraging over disorder, bilinear overlaps of these field act 
as order parameters at a saddle point. They are intro- 
duced at the cost of conjugated order parameters, that 
we indicate by a hat. The replicated free energy reads 
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Here \ a p = f3 2 J 2 p(p - l)q p af3 , a a/3 = a a 8 a p with a a = 
A aa (l - g Qa )/2 + l/(l - q aa ) and a a(i = p(3 2 J 2 (q a f i ) p ~ 1 . 
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Usually one assumes that no replica symmetry breaking 
occurs when one calculates the number of TAP states. 
However, since in eq. (^) it does occur at u = 1, it is not 
a well motivated assumption that it would not occur for 
u = 0. Having in mind the 1RSB solution of the ordinary 
partition sum, we are therefore led to investigate a simi- 
lar breaking scheme in the replicated TAP approach. We 
thus set q a/ 3 = {m a mp) = (q n - qi)S al 3 + qi£ a /3, where £ 
is Parisi's 1RSB matrix, having £ = 1 on all n/x diagonal 
blocks of size x*x, while zero elsewhere. A similar 1RSB 
pattern is assumed for the 5 other order parameters. At 
fixed x the 12-dimensional saddle point can be found ex- 
plicitly; details will be given elsewhere. As expected, the 
above replica expression for q is found back as solution 
of df(qo)/dqo = at qo = q. The result q\ — qo asserts 
that the mutual overlap between different states in the 
same cluster is equal to the selfoverlap. As in the replica 
calculation of the ordinary partition sum, x can still take 
any value. For the long time limit of the dynamical ap- 
proach the marginality condition should be taken, fl4| in 
the form given below eq. (||). It turns out that the fluc- 
tuation matrix then automatically has 3 zero eigenvalues, 
proving the marginality. The lowest acceptable value of 
F U (T; x) is obtained when x takes the boundary value 
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where one more eigenvalue vanishes. At the marginal 
point one then finds $„ = nI c /Nx(q, u), which implies 
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For u ~ 1 we thus recover the replica result (|7]) rather 
than the dynamical expression (^|). The latter expres- 
sion is still correct in the sense of CHS, namely as the 
component average free energy with entropy shifted over 
the constant 2c, so that it is continous at T c . However, 
after quantization T c then remains as residual entropy at 
T = 0. The derivation of eq. (|J) by CS is based on two 
assumptions that are both subject to criticism. 

From eq. (0) we observe that the complexity is tem- 
perature dependent: T(T) = I c /x(q(T)). Its increase 
towards smaller T shows that more and more compo- 
nents become available for the system. In other words, 
when decreasing the temperature phase space splits up 
more and more. At each temperature there appear new 
free energy barriers which can no longer be overcome. 
This well known notion of mean field spin glass theory, 
which nicely explains temperature cycling aging experi- 
ments, p^] is unknown in glass theory. Note that I{T) 
diverges as 1/T at low T. This leads to an infinite com- 
plexity at T — for time scales where the marginal state 
describes the physics. These will probably include the 
experimentally accessible time scales. 

It would be nice to have a measurable quantity that 
probes the multitude of solutions. One object that should 



be accessible, at least numerically, is the specific heat. 
The standard expression C = dll/dT = J2 a d(p a U a )/dT 
is likely to differ from the component average C = 
J2 a Pa,dU a /dT. [0 The interesting question is whether 
their difference is extensive. Based on experience in a 
toy model, we think it generally is. Since in the 
present model the energy fluctuations are too small at 
H = 0, Q it can only occur in a field. From the CHS 
or replica internal energy U in a small field we obtain 



l-C(T,H) = ±f3 2 J 2 (l + (p 
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On the other hand, a short calculation shows that C a — 
TdSa/dT remains only a function of q a at the marginal 
point, which in the present model is field-independent. 
This implies that C is field-independent as well, thus sat- 
isfying the Parisi- Toulouse hypothesis J|| . Interestingly 
enough, we find C < C, whereas Palmer derives the op- 
posite. Our reversed "dynamical" inequality is a new 
result that is due to the marginality. 

Some more aspects are of interest. Coming from the 
low temperature side, the transition has the nature of 
a continuous glassy transition. In a replica formulation 
states with breakpoint x = x(q) have massive longitu- 
dinal fluctuations (A' ~ T c — T) and massless transver- 
sal fluctuations. (The latter modes are commonly called 
the "replicon" ; we have proposed the more physical name 
"ergodon" jl2||). This structure is very reminiscent of the 
behavior of fluctuations in the SK model, be it much sim- 
pler. In an external field the degeneracy of the longitu- 
dinal eigenmode is lifted. It remains as a non-degenerate 
small eigenvalue which vanishes at the same point where 
the transversal "ergodon" becomes marginal. Therefore 
the same picture remains valid in an external field. In 
particular, there is a de Almeida-Thouless line where the 
paramagnet becomes unstable with respect to non-linear 
perturbations related to the 1RSB solution discussed. 

It would be useful to study the loop expansion. This 
may shed light on the lower critical dimension of the 
problem and on the relation with structural glasses. 

In finite (but large?) dimensions our explanation is 
incompatible with a droplet picture. Jl9[ The assump- 
tion that the condensed phase is non-degenerate implies 
T c /N = 0. Then it remains unclear from a thermody- 
namic point of view why a transition to a state with 
higher free energy can occur. It is also not obvious 
whether a droplet approach can explain our reversed in- 
equality C < C. 

In conclusion, we have considered a mean field model 
for a spin glass. As it undergoes a first order phase tran- 
sition to a state with one step of replica symmetry break- 
ing, it has many features of models for glassy transitions, 
possibly in the absence of disorder. Within this model 
we have found that the following scenario takes place: 
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- At T c the system freezes in a metastable state with free 
energy much higher than that of the paramagnet. 

- The component average entropy S drops at T c by an 
extensive amount X c , related to the degeneracy of the 
glassy states at T c . 

- The entropy S — S + X, the internal energy, and the 
free energy are continuous. 

- S vanishes at T = 0, so there is no residual entropy. 

- Below T c the glassy free energy lies below the param- 
agnetic one, and has a larger slope. Therefore the transi- 
tion is first order and unavoidable in the thermodynamic 
sense. The driving force is the complexity. 

- There is no latent heat, but a there is a drop in the heat 
capacity. 

- In an adiabatic cooling experiment the weights of the 
states are initially set at T c . However, there will be a 
further bifurcation of phase space at any lower T. 

- The inequality C/N < C/N is a direct result of the 
multitude of states and their marginality. 

Though having considered only one specific model, we 
are tempted to believe that the developed picture is very 
general. The role of disorder is probably not essential. 
What is needed is an extensive complexity. It seems to 
occur at many glassy and glass transitions. However, 
continuous transitions with infinite order replica symme- 
try breaking are known to behave differently. 
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